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PRECEDING  LAOS  hUU* 


AFSIM  —  An  Air  Force 
Satellite  Interactions  Model 


I  \  CYLINDER  MODKL  FOR  SCATHA 


1.1  Introduction 

Satellites  in  geosynchronous  orbit  at  approximately  G.  G  earth  radii  have  been 
found  to  charge  up  electrically  during  magnetic  substorms  (Deforest,  1  Rosen2). 
This  charging  process,  called  spacecraft  charging,  leads  at  times  to  discharges 
which  couple  radiated  energy  into  spacecraft  electronics,  leading  to  circuit  upsets, 
malfunction  of  satellites,  and,  on  occasion,  complete  cessation  of  satellite  opera¬ 
tion.  It  is  necessary  to  calculate  the  details  of  the  spacecraft  charging  process 
in  order  to  understand  the  physical  processes  taking  place,  to  design  spacecraft 
which  are  less  susceptible  to  charging,  and  for  the  purpose  of  taking  account  of 
potentials  in  the  satellite  sheath  in  the  analysis  of  particle  and  field  measurements 
by  instruments  on  board  satellites. 

The  first  theoretical  treatment  of  spacecraft  charging  at  geosynchronous  orbit 
was  carried  out  by  DeForest,  1  employing  a  current  balance  method.  The  satellite 
is  treated  similarly  to  a  Langmuir  probe  in  a  plasma.  Taking  into  account  back- 
scattering,  secondary  emission,  and  photoemission,  the  steady-state  floating 

(Received  for  publication  19  June  1979) 

1.  DeForest,  S.  E.  (1972)  Spacecraft  charging  at  synchronous  orbit,  J.  Geophys. 

Res.  77:651.  - 

2.  Rosen,  A.  (1975)  Spacecraft  Charging:  Environment  Induced  Anomalies. 

ALAA  Paper  75-91,  AIAA  Conference,  Pasadena,  CA. 
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potential  of  the  body  in  a  plasma  is  obtained.  In  this  method,  average  values  of 
secondary  emission  currents  and  backscattered  currents  are  employed. 

For  the  purpose  of  the  analysis  of  spacecraft  charging  and  design,  it  is 
necessary  to  be  able  to  compute  a  three-dimensional  distribution  of  potentials 
on  a  satellite  surface.  Inouye'*  developed  a  circuit  model,  in  which  currents  are 
injected  at  various  points  on  a  satellite  surface.  Approximate  expressions  for 
the  currents  are  obtained  from  probe  theory,  and  approximate  expressions  for 
secondary  emission  yields  are  employed. 

The  spherical  satellite  case  has  been  treated  in  the  steady  state  by  Parker.  * 

and  the  time -dependent  case,  by  llothwell  et  al'1  and  Hothwell  et  al."  The  method 

4 

of  Parker  consists  of  solution  to  the  l’oisson-Vlasov  equations  b\  the  following 
means:  a  series  of  iterative  operations  is  carried  out  in  which  ole.  troll  and  pro  on 
currents  to  the  spacecraft  surface  are  calculated  first;  then  after  the  deposited 
surface  charge  has  been  found,  the  Poisson  equation  is  solved  for  the  potential. 
Currents  to  each  surface  node  are  obtained  b\  trai  mg  particle  t  raioctorics  from 
the  node  into  the  surrounding  plasma  to  find  the  proportion  of  traiectories  which 
escape.  The  current  and  potential  computations  are  iterated  until  a  steady  state 
results. 

It  is  not  obvious  that  the  plasma  sheath  surrounding  a  satellite  has  a  steady 
state,  and  in  any  case  it  is  interesting  to  study  time-dependent  phenomena,  par¬ 
ticularly  the  emission  of  beams.  For  this  reason,  a  tune-dependent  spherically 
symmetric  description  of  satellite  charging  was  carried  out  b\  llothwell  et  al.  ’ 

The  method  of  computer  simulation  of  plasmas  is  employed,  til  which  the  motion 
of  thousands  of  electrons  and  ions  are  followed  in  their  self-consistent  electric 
field.  In  the  spherically  symmetric  model,  the  satellite  surface  is  considered  t. 
be  covered  w  ith  a  single  material.  Hackscattenng  and  secondary  emission  as  a 
function  of  incident  particle  energy  are  considered,  lti  addition,  photoemission  is 
accounted  for.  Active  control  of  the  satellite  potential  is  studied  bv  simulation  of 
the  emission  of  charged  particles  from  the  satellite  surface.  In  this  method,  an 

:i.  Inouye,  G.  T.  (1975)  Spacecraft  Potentials  m  a  Substorm  Fnvi  ronment ,  m 
Spacecraft  Charging  by  Magnetospheric  Plasmas,  A.  Hosen,  editor'. 

MIT  Press,  Cambridge,  MA. 

4.  Parker,  L.W.  (1975)  Computer  Method  for  Satellite  Plasma  Sheath  in  Steady- 

State  Spherical  Symmetry,  AFCH I .-TH -75-04 10,  AD  A015  066, 

5.  Hothwell,  P.  1..,  Ilubin,  A.  G. ,  Pavel,  A.  I..,  and  Katz,  I  .  (1975)  Simulation 

of  the  Plasma  Sheath  Surrounding  a  Charged  Spacecraft.  In  Spacecraft 
Charging  by  Magnetospheric  Plasmas,  A.  Hosen,  editor.  Mil  Press. 
Cambridge,  MA. 

ti.  Hothwell,  P.I..,  Rubin,  A.G.,  and  Yates,  G.  K.  (1977)  A  Simulation  Model  of 
rime -Dependent  Plasma -Spacer  raft  Interactions.  Proceedings  of  the 
Spacecraft  Charging  Technology  Conference,  C.  P.  Pike  and  It.  If.  l.ovell 
Editors,  AFGI.-TR-77-0051,  AD  A045  459. 


initial  configuration  is  sot  up  m  .1  spatial  grid  consisting  of  .1  sphcrit  al  satellite, 

.1  surrounding  thermal  plasma,  and  an  outer  boundary.  An  iterative  ca leulat ion 
is  carried  out,  m  which  each  of  the  plasma  particles  is  first  move.)  in  the  existing 
electric  field,  and  the  resulting  charge  configuration  is  employed  as  a  source  in 
the  calculation  of  the  electric  field  m  the  no\t  time  step.  The  computation  is 
stepped  forward  m  time  with  alternative  coin  put  at  ions  of  particle  motion  and  elec¬ 
tric  field. 

Kor  the  case  of  complex,  real  satellites,  it  is  desirable  to  take  into  account 
the  geometrical  details  and  materials  placements  and  properties  characteristic  of 
these  vehicles.  A  three-dimensional  quasi -statu-  code  called  NASP.AP,  develop.  ! 
In  Katz  ct  al,  accomplishes  this.  The  satellite  is  modeled  by  th<  1.  .  1. 

method  in  a  series  of  nested  grids,  each  of  which  has  a  grid  space  .  twice  as 
large  as  the  previous  grid.  Approximately  one  thousand  surface  cells  are 
accounted  for,  with  fifteen  materials  on  the  spacecraft.  Materials  properties 
such  as  backseattenng  and  secondary  emission  as  a  function  of  energy  as  accounted 
for,  as  is  photoemission.  The  external  magnetic  field  is  included  and  satellite 
spin  is  accounted  for.  The  photosheath  is  calculated,  as  well  as  charging  on  each 
satellite  surface  cell.  Trajectories  of  particles  from  emitters  on  the  satellite  are 
Computed  as  well. 

Because  NASt'Al’  does  not  treat  phenomena  01.  die  time  scale  of  plasma  oscil¬ 
lations  and  does  not  treat  space-charge  effects,  either  in  the  plasma  sheath  or  in 
emitted  beams,  we  have  developed  a  time -dependent  code  m  cylindrical  symmetry 
which  treats  these  effects.  This  code  is  designed  to  studv  beam -plasma  inter¬ 
action  effects,  active  control,  beam  neutralization,  sheath  stability,  space-charge 
effects  in  beam  emission,  and  satellite  surface  material  differential  charging. 

This  code  fills  the  gap  between  the  geometrical  simplicity  of  the  spherically 
symmetric  code  and  the  complexity  of  the  A  - 1  >  NASt'Al’  code.  It  has  the  capability 
of  handling  the  detailed  kinetics  of  spaoe-ehn rge  effects  omitted  from  the  NASt'Al’ 
code  and  geometrical  details  which  are  unresolved  in  the  SHEATll  spherically 
symmetric  code. 

1.2  \KSIM  Structure 

A  computer  program  has  been  developed  to  conduct  numerical  simulations  of 
plasma  interactions  w  ith  a  satellite,  assuming  the  satellite  to  be  represented  by  an 
infinitely  long  cylinder  and  treating  the  plasma  as  discrete  particles.  Vhc  system 
is  assumed  to  be  uniform  along  the  axis  of  the  satellite,  A  number  of  features  of 

7.  Katz,  I..  Harks,  O.  l\. ,  Wang,  S. ,  and  Wilson,  A.  (l!'77'  Oynannc  Modeling 
of  Spacecraft  in  a  Coll  is  toilless  Plasma.  Proceedings  of  the  Spacecraft 
Charging  Fechnology  Conference,  r.  P.  Pike  and  K.K.  1  ovell,  Editors, 

Al'C',1  -TK-77-OOT1I,  AH  A 04 a  4fd>. 
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'he  actual  satellite  arc  mi  -ot-po  rated  into  the  model  with  as  much  flexibility  as 
possible  to  ai  commodate  a  w  ide  range  of  conditions. 


1  t  triitlim  <>t  tin-  Model 

i.  t.  i  pa u  ru  i  i  Mo  floss 

Ihf  pa  ft  loles  a  if  1 routed  as  discrete  objects  lather  than  as  fluid  elements 
hut  ile  nut  necessarily  represent  individual  tuns  or  electrons.  The  current 
‘Install  associates  with  each  computer  particle  a  cluster  of  identical  electrons  01 
ions  so  that  the  several  million  particles  or  the  actual  plasma  are  represented  by 
a  few  thousand  computer  partii  les.  In  contrast  to  an  earlier  model  with  spherical 
symmetry,  the  parameter  relating  the  number  of  real  partii  les  for  each  computet 
particle  is  a  constant  foi  all  particles. 

The  particle  motions  are  determined  hi  the  electric  field  produced  l'\  the 
satellite  and  by  the  particles  themselves.  I' lie  program  allows  the  satellite  to  be 
charged  to  an  arbitrary  voltage  and  also  allows  the  potential  m  the  plasma  far 
from  the  satellite  to  be  set,  thus  influencing  the  rate  at  which  particles  emerge 
from  the  far  plasma  and  enter  the  region  near  the  satellite,  l'he  electi  n'  field  is 
calculated  from  the  electric  potential  obtained  by  solving  Poisson's  equation  m  a 
cy  lindrieal  geomet  ry  . 

I  hus,  the  equation  to  be  solved  for  the  potential  o  is: 


assuming  uniformity  is  the  /-direction.  With  a  coordinate  transformation 
S  In  r,  this  equation  becomes: 

2  2  2 

‘'1&  ,  a  f  *»'%> 

2  ■>  f 
US  dp  f 

which  is  the  form  for  rectangular  coordinates,  using  a  modified  density  r 

I  Ins  transformation  is  effected  by  choosing  the  radial  gridpoints  to  have  utii- 
fi'im  me  t  emetits  in  In  r,  rather  than  as  an  explicit  coordinate  t  rails  fop  mat  ion. 

A  standard  Poisson -solving  routine  (Hockney ‘S  for  rectangular  coordinate  systems 
is  then  employed  to  find  the  potential,  using  the  modified  charge  density,  l'lns 
approach  also  requires  an  auxiliary  routine  to  solve  for  the  potential  on  the  satel¬ 
lite  itself,  to  provide  a  necessary  boundary  condition  for  the  Poisson  solver,  l'he 

8.  Hockney,  H.W.  (1070)  l’he  potential  calculation  and  some  applications,  in 
Methods  in  Computational  Physics,  Vol.  0,  Academic  Press,  N.y. 


8 


auxiliary  routine  uses  a  moment  snu'S  expansion  to  del c rm me  the  satellite  po¬ 
tential  \n  terms  of  tile  charge  distribution  and  the  constant  value  assumed  for  the 
potential  m  the  far  plasma  region. 

No  magnetic  field  has  been  un  ,  rporated  into  the  program  at  this  tune,  hut  it 
would  he  feasible  to  add  a  magnet h  field  which  is  uniform  along  the  direction  of 
the  cylinder  axis. 

[’articles  are  lost'  from  the  system  by  striking  the  surface  of  the  satellite 
or  escaping  from  the  region  studied,  and  are  added  to  the  system  by  a  variety  of 
emission  processes  on  the  satellite,  described  below,  or  by  entering  into  the 
region  studied  from  the  surrounding  plasma. 

1 .  .1 .  2  SKCONOAin  I  MISSION 

Ions  or  electrons  which  strike  the  satellite  are  assumed  lo  generate  second. iry 
electrons,  taking  into  account  hoth  the  energy  and  angle  of  incidence  of  the  impact¬ 
ing  particle  to  determine  the  yield  of  secondary  electrons. 

J'be  equation  used  to  evaluate  the  secondary  emission  for  incident  electrons 
is: 


exp  If  1  1  -  cos  fl)  I 

w  here 

fi  number  of  secondary  electrons 
1'  energy  of  incident  electrons 

o  angle  of  incidence  with  respect  to  normal  to  surface. 

Vhc  parameters  l"  .  and  N  can  be  chosen  to  bo  different  over  the  angular 

o  1  o 

sectors  of  the  satellite  to  reflect  the  properties  of  different  materials  with  re¬ 
spect  to  secondary  emission. 

For  incident  tons,  the  equation  used  is: 

a  %  F.  see  o 

* 

w  here  a  (  and  Nj  reflect  properties  of  the  materials  in  each  angular  sector  of  the 
satellite. 

For  either  incident  ions  or  incident  electrons,  the  emitted  secondary  electron; 
arc  assumed  to  have  a  Maxwellian  veloeitx  distribution: 


S  l 


1  -  exp  (-IN  Fl") 


<\  I  > 


n-1 


The  velocity  dispersion  "b"  for  the  emitted  electrons  euti  depend  on  the  species  of 
incident  particle,  and  its  value  can  easily  he  changed  at  execution  time. 

The  direction  of  the  velocity  of  the  emitted  parti.  lev.  ha a  cosine  distribution 
in  angle  with  respect  to  tin-  normal, 

Ihe  corresponding  process  for  loti  emission  (sputtering!  has  not  been  treated, 
as  the  magnitude  of  that  process  is  relatively  small. 

1.  1.  d  HAt'KSt'AT  1’KKINC! 

Inelastic  collision  of  electrons  with  the  satellite  is  treated  as  a  hackscattering 
process,  again  taking  into  account  the  energi  and  angle  •  incidence  f  the  impact¬ 
ing  particle  to  determine  the  yield.  1'he  equation  used  here  is; 

.1  H  (A  K'nV'os  " 
o  o 

w  here 

number  of  hackscattered  electrons 
1!  energy  of  incident  electrons 

i>  angle  of  incidence  with  respect  to  normal  to  surface. 

I'he  parameters  A  and  11  can  likewise  he  set  to  different  values  on  the  different 
o  o 

angular  sectors  of  the  satellite. 

The  energy  of  the  emitted  particles  is  directlv  related  to  the  energy  of  the 

•(  ^ 

incident  particles,  according  to  f(v)  v  K“,  for  incident  energy  K,  provided 
•) 

m  v"  2  -  !•',  and  again  a  cosine  distribution  in  angle  is  assumed. 

As  with  the  sputtering  process,  hackseattcring  of  ions  is  neglected. 

1.3.4  IMIOTOTMISSU  \ 

The  emission  of  photoelectrons  from  the  satellite  surface  is  treated  in  a 
manner  similar  to  secondary  emission,  taking  account  of  the  incidence  angle  of 
solar  illumination  and  the  photoemissive  coefficients  of  the  different  sectors  of 
the  satellite  to  determine  the  photoemission  yield.  Thus,  the  relevant  formula  is; 

f  k  flcos  «») 
o 

w  here 

r  number  of  photoelectrons  per  unit  time  per  unit  surface  area 
f  solar  illumination  flux 

a  angle  of  solar  illumination  w  ith  respect  to  normal  to  surface. 


I  he  sola r  flux  is  essentially  zero  or  one,  depending  upon  whether  the  sector  is  in 

shadow  nr  riot,  while  k  is  a  photoeinissioti  coefficient  which  can  he  set  to  different 

o 

values  for  different  sectors  of  the  satellit. 

Currently,  the  velocity  spectrum  of  the  emitted  photoelect  runs  is  assumed  to 
he  a  Maxwellian  w  ith  the  cosine  angular  dependence  previously  assumed  for 
secondaries. 

1.  t.  ,  1‘AHTICl  i:  MKAMS 

<  or  espon ding  to  the  particle  beams  on  the  satellite,  there  is  a  routine  which 

emits  particles  from  a  localized  region  of  the  satellite  in  a  beam  pattern.  The 

beam  pattern  is  given  according  to  a  1  »  cos  "  "  law  ,  for  0  «  0  «  n  .  and  the 

O  o  o' 

veloc it \  distribution  of  the  beam  particles  can  he  chosen  to  he  either  monochro- 

2  2  2 

matic  (v  v  for  all  particles)  or  Maxwellian  < f < v )  v“  exp  (-v“  2v“)).  There 

o'  r  O 

are  effectivels  separate  particle  beams  for  ions  and  electrons,  and,  in  addition  to 
“  and  r  ,  the  location,  aiming  direction,  and  current  for  each  beam  can  he 
specified. 

An  option  is  available  by  which  the  total  charge  emitted  by  the  beams  can  be 
linked  to  the  total  satellite  charge,  thus  influencing  the  potential  on  the  satellite. 


2.  M.HTROS  I  \TI(  1*011  \  II  \l  <  M  i  l  I.VI  IONS  FOR  \  CVl  IMtRIC  Al. 

MOIIt  1 

2.1  Introduction 

Recent  investigations  of  spacecraft  interaction  with  the  surrounding  plasma 
have  led  to  the  development  of  a  two-dimensional  model  for  numerical  simulations. 
This  model  treats  the  spacecraft  as  an  infinitely  long  cylinder  and  studies  the 
motion  of  particles  in  the  cylindrical  annulus  between  the  spacecraft  and  an  arbi¬ 
trary  boundary  at  some  distance  into  the  plasma.  To  represent  the  motion  of  the 
particles  properly,  it  is  necessary  to  perform  calculations  of  the  electric  field 
which  is  established  by  the  distribution  of  charges  in  the  region  and  by  imposed 
potentials  on  the  spacecraft.  This  section  describes  the  potential  calculation 
method  that  has  been  implemented. 

2.2  Pumihi'i  Fipiatinn 

For  a  svstem  w  ith  translational  symmetry  in  the  Z-direction  (along  the  axis 
of  the  spacecraft  cylinder),  Potssor.'s  equation  is: 


l  iL  .  _L  '.f  £ 


(Ml 


(MKS  units) 


a here  the  electrostatic  potential  is  o,  and  p  is  the  charge  density  (per  unit  area). 
The  charges  are  effectively  infinitel>  long  rods  in  this  representation  and  interact 
via  a  logarithmic  potential. 

Converting  the  above  form  of  Poisson's  equation  to  suitable  form  for  a  dis¬ 
crete  grid  gives: 
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where  c  e(r  <•  ).  ,  p(r.  rt  1,  and  6  "  -0 

i.  J  »  J  ».  J  t  J  J  J-l 

I),  r  -  r(.  For  a  uniform  interval  in  the  radial  grid  points,  this  form  pro¬ 
duces  coefficients  of  ( (  ,,  which  depend  on  i.  However,  if  is  propor¬ 

tional  to  r  .  these  coefficients  are  .  oiistant. 
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Thus,  letting  .i  or.,  so  that  I) 
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o)r.,  Poisson's  equation  becomes: 
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This  formulation  corresponds  to  a  transformation  of  the  radial  coordinate  accord¬ 
ing  to  S  In  r,  which  gives 


^4  +  ££ 
as"  art" 


where  ,  r“  p  for  the  transformed  density  in  iS,  ■  1  coordinates. 

The  charge  density  is  letermined  from  the  number  of  charges  within  a  partii  - 

ular  grid  cell,  but  each  charge  is  treated  as  being  uniformly  distributed  over  a 

cell  in  (S.  »)  coordinates,  rather  than  (r,  H)  coordinates.  This  is  implemented  by 

2 

assuming  a  charge  density  distribution  ,  C  r"  for  each  cell,  normalized 

according  to  /  r  dr  drt  q  ,  and  then  using  the  average  of  r  p  over  the  cell  for 
J  t.t  i 

the  right  side  of  Poisson's  equation.  Thus,  r",(  i  is  replaced  by 

-q(  ^  (6  In  (1  +  u)l  in  the  formula  above. 

Note  that  different  charges  w  ithin  a  cell  do  not  interact,  and  only  the  net 
charge  within  a  cell  contributes  to  the  potential.  Thus,  the  grid  produces  a  poten¬ 
tial  that  is  both  "softer"  (reduced  short-range  interaction)  and  "smoother" 
(spatially  averaged)  than  the  actual  interparticle  interaction. 
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2.3  Round, ir\  Condition*  for  the  Model 


In  solving  Poisson's  equation,  the  boundary  conditions  must  be  specified  at 
some  radius  in  the  far  plasma  (Rp)  and  over  some  surface  on  or  within  the  space¬ 
craft.  The  spacecraft  itself  is  modeled  as  a  conducting  cylinder  completely  sur¬ 
rounded  by  a  dielectric  layer.  This  allows  the  charge  distribution  on  the  surface 
of  the  spacecraft  to  be  that  determined  only  by  particle  impacts  and  emission 
mechanisms,  yet  allows  for  the  specification  of  a  unique  fixed  potential  for  the 
spacecraft.  Conceptually,  therefore,  the  two  boundaries  are  concentric  cylinders 
with  a  constant  potential  on  each.  However,  neither  potential  is  necessarily  con¬ 
stant  in  time.  In  particular,  if  the  spacecraft  potential  is  not  fixed,  it  will  "float” 
with  respect  to  the  outer  boundary  potential,  depending  on  the  charge  distribution 
in  the  intervening  region. 

If  the  boundary  condition  is  left  solely  in  terms  of  the  potential  on  the  interior 
conducting  cylinder  (r  R.),  then  the  routine  for  solving  Poisson's  equation  must 
also  involve  the  conditions  at  two  interfaces:  the  conductor /dielectric  interface 
and  the  dielectric/vacuum  interface.  However,  if  the  interior  boundary  condition 
is  transformed  into  a  specification  of  the  potential  on  the  outer  surface  of  the 
spacecraft  (r  -  R^),  then  Poisson's  equation  need  only  be  solved  over  a  homo¬ 
geneous  region. 

This  transformation  can  be  accomplished  by  solving  for  the  surface  potential 
for  the  case  of  an  isolated  charge,  given  the  potential  at  the  outer  boundary  and 
the  specification  of  the  potential  on  the  interior  conductor  ("floating"  or  fixed 
potential).  The  solution  for  an  arbitrary  charge  distribution  is  then  the  appropriate 
superposition  of  solutions  for  isolated  charges.  (This  is  essentially  a  Green's 
function  approach. ) 

For  a  charge  q  at  (r  ,  6  ),  with  <j>  -  <t>Q  at  r  =  R^  (the  inner  conducting  cylinder) 
and  ij  at  r  Rp  (the  far  plasma),  the  potential  at  (R^,  0)  (on  the  surface  of  the 
spacecraft)  is  given  by: 
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In  principle,  the  summation  should  be  extended  to  n  op.  luit  on  a  discrete  grid  it 
in  l  ost  rioted  according  to  the  number  of  angular  sectors  !\  (  2 x  *1.  l  or  the  case 

of  a  floating  potential  on  the  spacecraft,  with  a  total  charge  Q  on  the  interior 
cylinder, 


♦„  0 1  ‘  :’«k, 


Some  modifications  to  the  above  expressions  are  required  to  accommodate 
the  convention  of  having  the  charges  distributed  over  each  cell,  rather  than  as 
discrete  charges.  Plus  is  accomplished  by  a  superposition  of  the  above  solutions 
for  a  charge  density 


In  ( 1  i  .>)  r 

over  the  region 
r 
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This  product's  a  potential  A(li  .  (I)  with: 
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where  a.  and  a.,  are  as  before. 
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The  condition  for  a  floating  spacecraft  potential  then  becomes 
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I'he  solution  is  now  completely  formulated,  given  a  method  oT  solving  Poisson's 

U 

equation  in  rectangular  eoontinates.  For  this  particular  application,  Hockney's 
rouricr  Analysis  Cyclic  Reduction  Method  was  employed.  I'he  Fourier  t ransform 
routine  for  Hockney's  programs  is  also  used  to  obtain  the  superposition  of  solutions 
for  setting  the  boundary  potentials  at  r  R  , 
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